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SUMMARY 

Prandtl's basic mixing length model was used to compute 22 of the 24 test cases 
for the Langley Working Conference on Free Turbulent Shear Flows. The calculations 
employed appropriate algebraic length scale equations and single values of mixing 
length constant for planar and axisymmetric flows, respectively. Good agreement with 
data was obtained except for flows, such as supersonic free shear layers, where large 
sustained density changes occur. The inability to predict the more gradual mixing in 
these flows is tentatively ascribed to the presence of a significant turbulence -induced 
transverse static-pressure gradient which is neglected in conventional solution pro- 
cedures. Some type of an equation for length scale development was found to be nec- 
essary for successful computation of highly nonsimilar flow regions such as jet or 
wake development from thick wall flows. 

SYMBOLS 

c p specific heat at constant pressure 



D nozzle diameter 

j geometry parameter, 0 or 1 for two-dimensional or axisymmetric flow, 

respectively 

k thermal conductivity 

l mixing length 

£ sd mixing length associated with second-derivative term in equation (5) 
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P 

Ap 

R 

Rb 


r 



T e 


mixing length constant 
Mach number 
Prandtl number, pcp/k 
Schmidt number 
static pressure 

static-pressure difference across shear layer 
constant in equation of state for a perfect gas 

Reynolds number based on characteristic body dimension, ii e u e r b lp e 
radial coordinate 

characteristic body dimension 
static temperature 


u mean velocity in streamwise direction 



u'v' turbulent shear stress 

v mean velocity in transverse direction 

v* fluctuating velocity in transverse direction 
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streamwise coordinate in physical plane 
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y transverse coordinate in physical plane 
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value of y where 



= 0.01, 0.05, 0.95, 0.99, 0.995, respectively 
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concentration 
mixing angle in wake flows 
ratio of specific heats 
width of mixing region 

eddy kinematic viscosity 
streamwise Von Mises coordinate 

momentum thickness 
dynamic viscosity 

density 

spreading parameter 
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u 2 

spreading parameter when — = 0 and M — 0 

U 1 

transverse Von Mises coordinate; or shear stress 
nondimensional transverse Von Mises coordinate, 

mixing angle for jet flows 


(p e u e ) 1/1+ir b 


Subscripts: 

c end of potential core 

<t center line 

e external stream 

i initial 

ip inflection point 

j jet exit plane 

l laminar 

n nth species 

p primary flow 

s secondary flow 

t turbulent 

1 high-velocity side of shear layer 

2 low-velocity side of shear layer 

I, II regions I and II in jet flows (figs. 2 and 4) and wake flows (fig. 3) 

TIT region III in jet flows (figs. 2 and 4) 
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Superscript: 


j geometry parameter where j = 0 in two-dimensional flows and j = 1 in 

axisymmetric flows 

THE MIXING LENGTH CONCEPT 

The classic phenomenological model for the turbulent shear stress is the mixing 
length theory (refs. 1 and 2) developed by L. Prandtl in 1925. Although the model lacks a 
rigorous physical basis, it has nonetheless proved to be quite successful in both boundary- 
layer and free shear flow calculations. 

In formulating his model, Prandtl assumed that the Reynolds stress Is produced by 
momentum transfer normal to the mean flow direction from regions of high momentum to 
regions of low momentum. The mixing length l is then defined (in rough analogy to the 
mean free path of a molecule) to be the average distance over which a fluid element trans- 
fers momentum. A detailed description of the model’s theoretical development is avail- 
able in many texts (e.g., ref. 3, p. 277 and ref. 4, p. 545). 

Mathematically, for quasi -parallel shear flows the eddy viscosity is expressed as 
the product of the square of the mixing length and the absolute value of the local mean 
velocity gradient, that is, 


e = l 


2 


9u 

9y 


(i) 


Therefore, the turbulent shear stress for two-dimensional or axisymmetric flow becomes 


r t = -p u V = pe — = pi 2 

9 y 


9u 

9y 


9u 
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( 2 ) 


where the mixing length is generally taken to be the product of an empirically determined 
constant and some characteristic scale of the mixing region, that is, 


l = Constant x 6 


or 


Constant = — 
6 


( 3 ) 


71 


Thus, 


/if 62 
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9y 


au 
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( 4 ) 


Equation (4) has been widely used in the calculation of a variety of flows with both 
analytical and numerical solution methods. After early use in pipe flow and flat-plate 
boundary -layer analytical solutions, the mixing length model has been particularly suc- 
cessful in the numerical computations of both incompressible (e.g., refs. 5 and 6) and 
compressible (e.g., refs. 7 and 8) turbulent boundary -layer flows with appropriate defini- 
tions of the mixing length specified for the "wall" and "wake" regions. For example, 
Patankar and Spalding (ref. 7) defined the mixing length as a numerical function of the 
distance from the wall and the boundary-layer thickness. Bushnell and Beckwith (ref. 9) 
then extended the use of the mixing length model to include the compressible nonequilib- 
rium case. 


A number of well-known approximate analytic solutions were developed for low- 
speed free flows by using Prandtl's mixing length model. Prandtl studied the smoothing 
of a velocity discontinuity (ref. 4, p. 687) and assumed similar velocity profiles. Tollmien 
assumed similarity and presented solutions for both the circular jet (ref. 4, p. 699) and the 
two-dimensional jet (ref. 4, p. 696) issuing into still surroundings as well as the two- 
dimensional two-stream mixing layer (ref. 4, p. 689). Schlichting (ref. 4, p. 691) in his 
thesis first examined the two-dimensional wake. Using measurements in the wake of a 
circular cylinder, he found 1 / 6 to be 0.18. Kuethe (ref. 10) extended Tollmien's solu- 
tion for the mixing region of a plane jet to the general case of the mixing of two parallel 
streams with different velocities. Kuethe also solved the problem of the mixing region 
surrounding an axisymmetric jet issuing into still surroundings by assuming local simi- 
larity. Using experimental measurements and taking the mixing length to be proportional 
to the width of the mixing region, he obtained 1/6 = 0.0705. Squire and Trouncer (ref. 11) 
considered axisymmetric coaxial jets and used l/6 = 0.082. 

Several objections to the mixing length concept have been raised. For instance, 
according to equation (2) the eddy kinematic viscosity € always vanishes at points of 
zero velocity gradient, a condition which is physically unrealistic (e.g., the center line of 
a jet or wake). Prandtl (ref. 1) recognized this difficulty and proposed that e be con- 
sidered proportional to a statistical mean of 9u/9y in the vicinity of a velocity maximum. 
The mixing length model for shear stress thus becomes 


.2 9u 

T t = ^ * 



ni/2 


(5) 
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where Z g{ j is an additional length which must be determined experimentally. Equa- 
tion (5) has seldom been used in analytical solutions because of the mathematical diffi- 
culties of including its more complex form. (There is no such problem if numerical 
solutions are employed.) 

Another deficiency of the model given by equation (4) becomes evident in the study 
of the two-stream mixing layers. In this case the analytic solution exhibits a discontinu- 
ity in a^u/Sy 2 at the mixing region boundaries instead of predicting the expected asymp- 
totic approach of the mixing region velocity to the two free-stream values. Equation (5) 
is free of this problem. 

In addition, the original assumption of constant mixing length across the mixing 
region has led to inaccurate predictions of the fluctuating components, according to 
Liepmann and Laufer (ref. 12), despite generally good agreement of the mean velocity 
profiles with experimental data. 

Later Alexander, Baron, and Comings (ref. 13) demonstrated that Prandtl's original 
development required the assumption that the turbulence intensity and mixing length be 
isotropic. They concluded, however, that this first-order assumption did not seriously 
limit the application of the model. The original physical concept of a mixing length seem- 
ingly requires that its magnitude be small in comparison with mean flow dimensions. 
However, the fact that it has been found experimentally (ref. 14) that the mixing length is 
not small in comparison with mean flow dimensions is often considered a basic defect of 
the mixing length theory. 

Recently several authors (e.g., ref. 15) have shown, by means of the turbulent kinetic 
energy approach, that the mixing length hypothesis implies that production equals dissipa- 
tion of turbulence energy. Comparison of figures 6-9 and 7-44 in Hinze (ref. 3) shows 
that this assumption is much more defensible in boundary layers than in free shear flows. 

Most investigators abandoned the mixing length in favor of Prandtl's later constant 
exchange coefficient eddy viscosity model (ref. 16), whose mathematical definition was 
more easily incorporated into analytical solutions. This "new" model overcame some of 
the difficulties of the mixing length model and gave slightly better agreement with experi- 
mental data. It has become historically the most widely used incompressible model. 

Recently Elassar and Pandolfini (ref. 17) numerically investigated the compressible 
free shear layer in a backstep region and concluded that 1/6 must have different values 
in the similar and nonsimilar regions. Wagner (ref. 18) used the mixing length model with 

the second-derivative term (eq. (5)) and I = ^ = 0.08 to correctly predict the develop- 
ment of the hypersonic turbulent wake produced by a wedge. 

The first attempt to systematically apply the mixing length model to a significant 
range of free turbulent flows was Harsha's study (ref. 19) of eddy viscosity models. 
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Harsha applied a single mixing length constant, 0.082 (the value used by Squire and 
Trouncer (ref. 11)), to both compressible and incompressible flows. Apparently, for all 
regions except the potential core of a jet, he defined the width parameter to be twice the 

half-velocity radius ^the radius for which u = ^u^ + Ugj^j. In the shear-layer region of 

jets he used the width of the mixing region. This approach led to poor overall agreement 
with the selected experimental data although he obtained good prediction of the velocity 
potential core length for the subsonic circular jet and the subsonic coaxial air-air mixing 
cases. His predicted potential core was too short for the compressible jet issuing into 
still air and too long for hydrogen-air coaxial mixing flows with mass flow ratio P e u e /pjUj 
greater than 2.0. For all these calculations, the asymptotic slope and the center-line 
velocity decay curve was underpredicted. However, the mixing length model did predict 
the early center-line velocity increase for both an axisymmetric and a two-dimensional 
subsonic wake very satisfactorily, although the downstream values were underpredicted 
in both cases. 


THE "RATIONAL” APPROACH 

The present effort was the result of a desire on the part of the conference organiz- 
ers to have the mixing length eddy viscosity model compared with newer formulations. 
Recent success (refs. 9, 20, and 21) at Langley in the use of Prandtl's model in high- 
speed boundary layers indicated that the model should be given further consideration in 
free shear flows as well. Rather than repeat the Harsha approach (ref. 19) with the 24 
conference test cases, an approach which might be termed a more "rational" approach 
was selected in which each class of flows would be given individual consideration. The 
designation "rational" merely indicates that the approach was designed to better repre- 
sent the actual physics of free mixing than can be accomplished by conventional use of 
the basic model. 

On the basis of boundary-layer experience, the so-called "rational" approach was, 
at the beginning of the investigation, visualized as the inclusion of the following modifi- 
cations in the basic mixing length model (eq. (2)): 

(1) A consistent set of values for the proportionality constant If 6 were to be used. 
Originally it seemed that this might involve a different constant for each region of a jet 
or wake. However, the results to be presented in this paper indicate a single value of 
1/6 for all two-dimensional flows and a single value for all axisymmetric flows, a con- 
clusion in agreement with the experience of Spalding's group at Imperial College (refs. 15 
and 22 and paper 11 of this compilation). 


74 


(2) Most importantly, a simple length scale equation would be included to provide a 
reasonable transition from the relatively small scale of turbulent motion in the shear 
layer of a jet to the larger scales of turbulence in the far field. Similar consideration 
would be given to wake and coaxial flows developing from wall flows. This approach had 
been successfully used in the slot injection studies of references 20 and 21. 

(3) Recent experience (ref. 23) in turbulent-boundary -layer calculations has shown 
a pronounced increase in turbulent shear stress near the end of transition and the begin- 
ning of turbulent flow. This "low Reynolds number effect" was to be investigated for any 
relevant free shear flow test cases. However, for the test cases considered herein this 
effect was not included due to the lack of information concerning possible magnitudes and 
variations of the "low Reynolds number effect" in free turbulent flows. 

(4) The use of an intermittency factor in the calculation procedure was to be 
explored. Wagner (ref. 18) had previously used intermittency in his hypersonic wake 
calculations. 

(5) The effect of including the second-derivative term in the mixing length model 
(i.e., eq. (5)) would be studied for some of the test cases. 

SOLUTION PROCEDURE AND TURBULENCE MODELS 
Basic Equations 

Solutions were obtained by using a computer code developed by Sinha, Fox, and 
Weinberger (refs. 24 and 25) for chemically nonreacting, quasi-parallel shear flows. 

The boundary -layer equations with suitable boundary conditions are assumed to 
describe the motion of free shear flows. Only quasi-parallel flows are considered. 
Effects of longitudinal curvature and transverse pressure gradient are not included. For 
the mixing of dissimilar ideal gases the governing equations are as follows: 

Conservation of mass 

^( puyi ) + |;( pvyi ) = 0 <6a) 


Conservation of streamwise momentum 


pu 


0X By 9x yj 9y' v ByJ 


(6b) 
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Conservation of energy 


9T 0p e 1 0 L i 0T 


9y 9x y j 0y\ 9y 


ii-fkyi *1\ + mH (c„ n — — 

..i 9vl 9y/ \9y/ N c „ ^ | p ’ n 0y 0y y 


Conservation of nth species 


0a n 9a Ida i 8 «n 
pu —R + pv — = - — yJ — 


0 y v j 8 y\ N sc 0 y 


where j = 0 or 1 for two-dimensional or axisymmetric flow, respectively. The associ- 
ated boundary conditions are as follows: 

At the center line (y = 0), 

. *1 . !2S . o <6e> 

9y 0y ey 

for all values of x. 

At the edge of the mixing region, 


lim u(x,y) = u £ 


lim T(x,y) = T e 
y-°° 


lim a n (x,y) = a n e 
y— 00 


Equations (6a) to (6f) are solved in the modified Von Mises plane defined by the fol- 
lowing transformation: 

(1 + j)^ = Puyj (7a) 

9y 


(1 + j) T] = -Pvy J 
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and 



In the transformed plane only the momentum, energy, and species equations are 
needed since the continuity equation is inherently satisfied in the definition of the trans- 
formation. The transformed equations are then nondimensionalized with respect to 
external stream conditions and r b , a characteristic length. It should be noted that 
Wagner (ref. 18) has shown the analysis of references 24 and 25 to be incorrect for longi- 
tudinal pressure gradient flows since the continuity equation was not satisfied by the non- 
dimensional Von Mises transformation used in that analysis. The correctly transformed 
equations of motion are as follows: 

Momentum 



Energy 


9T 


J-li . M i?( lnu e) n + 4! 4/ la. a iOl cr *l) 

l+j\ / d£ 3r c p 9 ?\N Pr f j P QtJ 


+ 





(9b) 


(9c) 


Equations (9a) to (9c) were applied to turbulent flow by assuming that 


M = H + Mt 


( 10 ) 
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and by using turbulent Schmidt and Prandtl numbers. Suitable expressions were used for 
the laminar viscosity of the various gases used. 

The flow was computed by using an implicit finite -difference technique (which is 
described in detail in refs. 24 and 25) to march downstream from the initial input profile. 

Use of Von Mises Transformation 

Apparently, most of the available computer codes for quasi -parallel free viscous 
flow solutions (e.g., refs. 7, 24, and 26) solve the boundary -layer equations in a Von Mises 
plane rather than the physical plane. For programing simplicity, a constant step size in 
the direction normal to the flow is usually assumed as in the procedure of reference 24. 
From equation (7a) this radial coordinate t in the method of reference 24 is defined by 



Equation (11) shows that specifying equal At steps implies specifying equal increments 
of axial-direction mass flow across the computed region. Thus, larger steps (in terms 
of the physical plane) are located in regions of relatively low velocity (with respect to the 
free-stream or jet center-line velocity); this significantly reduces the accuracy in the 
crucial mixing regions. This problem is particularly evident in the calculation of com- 
pressible flows, especially the shear-layer calculations. As an example, in a typical 
Mach 5 two-dimensional shear-layer calculation with \i2/ u l = 0-05, 49 of the 90 points 
across the mixing region profile are located in the region where u is approximately 
constant (u/uj = 0.90) . 


Turbulence Models 

In accordance with the "rational approach" previously outlined, a conceptually sim- 
ple (i.e., first order) model for the length scale variation in the streamwise direction was 
formulated for each of the major classes of flows. Thermal transport and mass transport 
were related to the momentum transport through the use of constant values of the turbulent 
Prandtl number and Schmidt number, respectively. 

Two-dimensional shear layers . - Figure 1 shows the initial velocity profile used to 
generate the two-dimensional shear layers of test cases 1, 2, and 3. The same type of 
step profile was used for both the initial temperature and initial concentration profiles. 
Because of the inaccuracies introduced by the Von Mises transformation at the edge of the 
mixing region, the following definition for the width 6 of the region was necessary: 

6 = (Scaling factor)^ g5 - y ^ (12) 
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where y g5 indicates the point in the mixing layer where - u^j^u e - u<g = 0.95. The 
scaling factor is necessary to obtain the "true" value of 6 (i.e., y gg - y and was 

taken to be 1.425 from the mean fit curve to the experimental data for two-stream mixing 
given by Halleen (ref. 27). This definition of width was used also for test case 5 and in 
the potential core region of all jet flows. 

Since the initial boundary-layer profiles were given for test case 4, the two- 
dimensional version of the "coaxial jet length scale model" was applied to that case. 

Jets into still air .- Figure 2 shows the length scale model developed for jets issuing 
into still air. The nondimensional mixing region width 5 is partially defined by the mix- 
ing angle <fi. This definition allows a smooth transition between the relatively small mix- 
ing scale at the end of the potential core (which is designated region I) and the large scale 
turbulence downstream (region III). 

For region I, Bj is given by equation (12); for region II, 

5 n = (Scaling factor)^ g5 - y^ + ^x - x c jtan $ (13) 

where x c is the location of the end of the potential core. Region II extends to the point 
downstream at which 6 jj becomes 2(Scaling factor)^ g ^ - y^. Finally, in region in, 

B m = 2 (Scaling factor)^ g& - y t ) (14) 

The scaling factor for regions II and HI was taken to be 1.2 from the similarity profile for 
a self-preserving jet given by Wygnanski and Fiedler (ref. 28). 

Since the governing equations were nondimensionalized with respect to the external 
streams, some small arbitrary edge velocity was necessary. However, the velocity ratio 
u^u e at the initial profile had to be large enough to reasonably represent the physical 
problem being studied but not so large that the accuracy on the low-velocity site would be 
limited by the peculiarities of the Von Mises transformation discussed previously. The 
values of u^u e used in the initial profiles for the jets issuing into still air are listed in 
the following table: 


Test case 

Initial u^u e 

6 

34.65 

7 

35.37 

8 

41.80 

18 

70.0 

19 

55.30 


79 



This length scale model was also used to compute test case 13, a two-dimensional 
jet in a moving stream. 

Wakes .- Figure 3 shows the length scale model employed in wake flows. Again, a 
linear variation defined by the angle 0 is assumed to approximate the increased scale 
of turbulent motion when the boundary layers on the wake -generating body mix in the ini- 
tial wake region. For region I in figure 3, the width 6 is given by 

6j = (Scaling factor) (y.95-y<t) + ( i - i H) tan ' 3 (15) 

where x^ is the location of the start of region I which is usually the trailing edge of the 
body. Region I terminates when the two terms in equation (15) are of equal magnitude. 
The width 6 in region II is then given by 


6ji = 2(Scaling factor)^ .95 - y*) (16) 

The same scaling factor, 1.2, utilized in the jet flows downstream of the potential core 
was also used for wake flows. 

Coaxial jets .- The model shown in figure 4(a) was employed for coaxial jet flows. 
The primary jet flow is assumed to mix with the initial boundary -layer flow in a region 
whose growth is characterized by mixing angles <p (p 2 , and <p g. In all calculations 
these angles were assumed to be equal. This model results in a radial distribution of 6, 
such as the one shown in figure 4(b), and features two overlapping scales of turbulent 
motion, 6 0uter and S Inner , given by 


6 t = 6 . 
Inner p,i 


2 x 


-Xi) 


tan cp 


(17) 


6 Outer " (y.995 “ r j 




tan 4> 


(18) 


where x^ is measured from the exit plane of the primary jet nozzle. In the calculations, 
6 is 6 Inner for y < r p and S 0uter for y > r g . Between r p and r s , 6 is found 

by a linear interpolation between 5 Inner and 6 0uter . The nondimensional radii ? p 
and r s are defined by 


r P"2 


r i + 


^x - x^tan <fi 


(19) 
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and 




r s = 


f j * | 5 S , i 

|l[ f i + v + ( 5 


x - xj] tan ( p 


for x = x r 


for x > x f 


( 20 ) 


where ?j is the inner nozzle radius, x c is the location of the end of the core region, 
and 6 S j is the initial thickness of the external boundary layer. This length scale 
modeling allows (to zeroth order) for the orderly growth of the initial boundary layer (or 
shear layer) scales. Far downstream the more conventional approach of a single large 
scale for the entire flow is recovered. 


A differential equation for the local length scale would obviously be better than the 
present approach just outlined. The purpose here is simply to point out that if the physics 
of scale development is taken into account, predictions are considerably improved. In the 
present approach, this scale development is handled in the simplest manner possible. 


POSSIBLE INFLUENCE OF TRANSVERSE STATIC -PRESSURE 
GRADIENT ON TURBULENT FREE MIXING 


The conventional approach to turbulent free mixing calculations (the approach used 
in the present paper) involves solving the boundary-layer equations with some turbulence 
closure model. On the basis of an order of magnitude analysis, the normal (or trans- 
verse) static-pressure gradient is assumed to be quite small in the basic boundary -layer 
equations. By assuming constant static pressure in the transverse direction, the solution 
of the normal momentum equation is unnecessary. Recent studies (refs. 29 and 30) indi- 
cate that this assumption of constant static pressure is not correct in boundary layers 
where there are large density changes across the layer. Changes in static pressure of 
10 to 30 percent have been observed in nominally two-dimensional, high Mach number 
(M > 6) boundary layers. 

The normal momentum equation for two-dimensional turbulent quasi -parallel shear 
flows takes the form 
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where density fluctuations have been neglected in an order of magnitude analysis. Hinze 
(ref. 3) notes that for low-speed flows the normal static -pressure changes induced by the 
turbulence (from eq. (21)) were 2 orders of magnitude larger than the static-pressure 
changes induced by the mean flow. The basic question is whether or not the normal 
static-pressure gradient induced by the turbulence (eq. (21)) can become sufficiently 
large to affect the spreading rate and turbulence structure of a free turbulent mixing 
flow. The success of various calculation methods (turbulence closure models) in low- 
speed flows indicates that the influence is most likely negligible in these flows. However, 
large density changes would presumably increase the magnitude of the turbulence-induced 
normal static -pressure gradient (from eq. (21)); therefore, the answer to the basic ques- 
tion may well be that the spreading rate and turbulence structure are in fact affected 
where large density changes occur. 

Estimates were made of the probable static-pressure change for an Mj = 5 free 
shear layer obtained recently in the nozzle test apparatus at the Langley Research Center. 
(These data are discussed in paper 2 by Birch and Eggers.) By using the integrated form 
of equation (21) and estimates of v'2 from the subsonic shear layer of Liepmann and 
Laufer (ref. 12), the change in static pressure for this flow would be of the order of a 
100-percent decrease near the maximum shear region. Obviously v'2 and perhaps the 
Reynolds stress must be lower than the normalized values of reference 12 for this Mj = 5 
case, since such a transverse static-pressure gradient would not be tolerated by the flow. 

There are at least two possible mechanisms connected with a turbulence -induced 
3p/9y which could affect the spreading rate of a free turbulent shear flow: 

(1) Direct influence of turbulence -induced 9p/9y upon the normal mean velocity v 
and hence directly upon the spreading rate (refs. 31 and 32). A pressure gradient dp/dx 
is imposed within the shear flow. 

(2) Indirect influence of turbulence -induced 9p/9y upon v'2 and u'v' by means 
of a turbulence field adjustment to decrease 9p/9y until a balance is achieved. 

Obviously both of these mechanisms as well as others could be operating simultaneously. 
The important point is that as the density change across a turbulent free mixing region 
increases, the influence of the turbulence -induced 9p/9y may become increasingly more 
important, especially in determining the absolute physical dimensions of the flow field 
(i.e., the entrainment rate). 

Transverse static-pressure data for free shear layers are available in references 32 
to 34. For low-speed flows the data indicate a minimum pressure of the order of 2 per- 
cent near the inflection point in the velocity profile. The low-speed data of Lee (ref. 35) 
and Jones and Spencer (ref. 36) indicate small static -pressure changes (less than 1 per- 
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cent) for shear layers with 


u 2 

— = 0.35 and 0.30, respectively. 


Further data for the trans- 


verse static-pressure distribution in jets are available in references 31 and 37 and some 
data for wakes are found in reference 38. The available data do seem to indicate an 
increase of transverse static-pressure gradient with increasing Mach number across the 
shear layer. 


It is interesting to note that the integrated form of equation (21) is 


or 


Ap 



Ap 

P t 


v' 2 y <t M 2 R( t T <t 

u t 2 >ip ^ ^p T ip 


( 22 ) 


where the subscript ip indicates the inflection point in the shear-layer profile. There- 
fore, Mach number may well have an effect greater than simple density changes (caused 
by gas composition) upon possible changes in entrainment rate due to 8p/9y. 

The possibility that static-pressure variation may be responsible for the rather 
large changes in spreading rate with Mach number noted in connection with test case 2 in 
the present paper is entirely speculative at this point. There are several possible 
approaches that could be used to investigate this possibility. One such method would be 
to incorporate the complete normal momentum equation in a m etho d such as Donaldson's 
(ref. 39) where v' 2 is directly computed. If an equation for v'2 is lacking, it could be 
related to the Reynolds stress by some constant or function. Again, it would be necessary 
to solve the normal momentum equation (including the convective terms). The best method 
(i.e., most correct) of including the normal momentum equation would be a full field 
Navier-Stokes solution with turbulence terms included. However, if molecular diffusion 
of momentum is neglected in the normal momentum equation, a first-order solution may 
be obtained numerically by including the equation in a streamwise iteration loop with the 
usual boundary -layer equations. 

There are, of course, other mechanisms which could also account for the disagree- 
ment between experiment and theory for case 2. These possibilities include the increased 
influence of p'v 1 and other terms in the second-order correlation equations. These 


/ 
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terms are usually negligible at low speeds but presumably become more important at 
high speeds. Currently there are no available models or data for these terms in high- 
speed flows which could be used to check this hypothesis. 

RESULTS AND DISCUSSION 

Solutions were obtained for all 17 of the primary test cases and for five of the seven 
optional cases. The two cases omitted were case 22, a hydrogen-air coaxial jet flow, and 
case 23, a subsonic compound coaxial air jet. These cases were omitted because of 
numerical and programing difficulties encountered rather than failure of the basic turbu- 
lence models being used. 

During the investigation it was found that 1/6 could be considered a function only 
of flow geometry. A value of 0.07 was used in all two-dimensional flows and 0.05 was 
used in all axisymmetric flows including the shear-layer region of jets. 

In general, for constant density subsonic flows the turbulent Prandtl and Schmidt 
numbers were taken to be 1.0. For variable-density flows either 0.7 or 0.8 was used, 
with the Prandtl number assumed equal to the Schmidt number in all cases. 

Two-Dimensional Shear Layers 

Test case 1.- Two-dimensional subsonic constant-density free shear layers were 

computed for velocity ratios U 2 /uj of 0.05, 0.2, 0.4, 0.6, and 0.8 by using "step" input 

profiles for velocity and temperature as previously described. Since the equations were 

normalized with respect to the secondary stream, a nominally zero value of U2 was not 

possible in the present program; — = 0.05 was used in both test cases 1 and 2. The 

U 1 

results from the computations of test case 1 are shown in figure 5. The spreading 
parameter a is defined by the following relationship: 



u - u 2 

where Ya and Yd are the lateral distances between points at which is 0.1 

A iD Uj - U2 

and 0.9 at longitudinal stations X A and X B , respectively. The spreading parameter o 
is thus the reciprocal of the spreading rate of the shear layer. 

Equation (12) was used to define the width 6, and 1/6 was taken to be 0.07. By 
extrapolating the computed values of a to the case where U 2 /uj is zero, o 0 was 
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found to be 11.4. This value was used to normalize the computed values of a. As shown 
in figure 6, the results agree with the well-known relationship (from ref. 40) 


\ 


% u - u 2 
° *1 + u 2 


Test case 2 .- Two-dimensional shear layers were computed for five values of Mj 
to determine the effect of Mach number on the spreading rate (test case 2). The results 
are shown in figure 6 along with experimental data from references 12, 28, 31, 34, and 41 
to 50 as well as the unpublished data recently obtained at the Langley Research Center at 
Mi = 5.0. These shear-layer data have been adjudged the most reliable by Birch and 
Eggers in paper 2 of this compilation. 

Equation (12) was used to define the width 6. With an 1 / 6 of 0.07, only a slight 
increase in o is predicted for increasing values of Mi in striking contrast to much of 
the available data. At Mi = 5.0, the effect of decreasing the mixing rate by lowering 1/6 
is shown. The cause underlying the inability of the mixing length model to predict the 
spreading rate of supersonic shear layers is perhaps directly related to the validity of 
the fundamental assumptions inherent in the boundary -layer approach to free mixing prob- 
lems. One possible explanation is the existence of a transverse static-pressure gradient 
which could affect the turbulence spreading rate. The effect of such a pressure gradient 
was previously discussed in detail in this paper. There is also the possible influence of 
"low Reynolds number effects" as discussed by Birch and Eggers in paper 2. The data 
that may be influenced by such effects are indicated in figure 6. It should also be noted 
that the values of a computed with ^ = 0.05 were scaled to correspond to ^ = 0 

(using fig. 5). 


Test case 3.- Test case 3 was designed to determine the effect of density ratio on 

° Un 

the growth of fully developed two-dimensional shear layers. With = 0.2, calculations 

U 1 

were made for shear layers with density ratio P\jP2 approximately 1/14, 2, 7, and 14. 

Pi ! ' 

For all but the — = — shear layer, the density difference was accounted for in the pro- 
P 2 14 

gram by two methods: the subsonic mixing of gases with different molecular weights, and 

the mixing of two air streams with supersonic or hypersonic velocity in the primary 
Pi 1 

stream. The — = — shear layer was obtained by mixing hydrogen in the primary 
P% 14 

Pi 

stream with air in the secondarv stream. The — = 2, 7, and 14 Shear layers were the 

p 2 

result of mixing hydrogen and helium, helium and air, and air and hydrogen; the corre- 
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sponding air into air equivalent density ratios were computed by using = 2.24, 5.48, 
and 8.06, respectively. Turbulent Prandtl numbers of 0.7 and 0.8 were used for the dis- 
similar gas and supersonic flows, respectively; the turbulent Schmidt numbers were also 
0.7 and 0.8, respectively. Appropriate input profiles of the type shown in figure 1 were 
used in each instance. The results, shown in figure 7, indicate that a is a similar func- 
tion of density differences irrespective of the mechanism producing that difference. How- 
ever, it should be remembered that in test case 2 the supersonic results were greatly in 
error when compared with the data. 

Test case 4 .- Test case 4 considers the subsonic mixing of two streams initially 
developed on a symmetric airfoil with a 10° trailing edge. The coaxial jet length scale 
model was used for this calculation. A representative small angle of 8° (the value used 
in the slot injection studies of refs. 20 and 21) was used. The predicted profiles at two 
locations, x = 12.7 cm and x = 76.2 cm, are compared with Lee's experimental mea- 
surements (ref. 35) in figures 8(a) and 8(b), respectively. For this test case and test 
case 5, the computed profiles are matched with the experimental profiles at the point 
where u = ^uj_ - u 2j ^ 2 since the numerical calculation does not satisfy the proper bound- 
ary conditions (ref. 51). The predicted profile at x = 12.7 cm indicates only slightly 
more mixing than the experimental data. Much farther downstream at 76.2 cm the pre- 
diction indicates significantly more mixing. Better agreement could have been achieved 
by adjusting the angles 0j, 02 » and 03 but such manipulation was beyond the intent of 

the current investigation. It should be noted that the streams mix with a fairly large ini- 
tial angle (10°). 

Test case 5 .- The Hill and Page (ref. 50) supersonic = 2.09j shear layer (test 
case 5) was computed with 1/6 = 0.07, where 6 was defined by equation (12). Computed 
profiles at two stations are shown in figure 9. The effect of the different boundary condi- 
tions used in the computations is evident on the low-velocity side of the shear layer. The 
data were obtained in a cavity-type flow and may be influenced by the "low Reynolds num- 
ber effects." Both calculated profiles indicate more rapid mixing than the experimental 
data, which is the expected behavior from the results of test case 2. (See section entitled 
"Possible Influence of Transverse Static -Pressure Gradient on Turbulent Free Mixing.") 

Jets Issuing Into Still Air 

For all flows of this type the length scale model of figure 2 was employed. An 1/6 
of 0.05 was used in all three regions of the axisymmetric jets including region I. The 
value of 1/6 should actually be varied from 0.07 at the nozzle exit where the shear layer 
is thin and essentially two-dimensional to 0.05 at the end of the core where the flow is 
fully axisymmetric. Hence, for flows with relatively short potential cores (e.g., test 
cases 6 and 8) the use of a single value lower than the average value through the core 
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results in less mixing and hence too long a potential case. The present model predicted 
the decay of the center-line velocity remarkably well by using an angle <fi of 8°. 

Test case 6.- Figure 10 shows the results of the computation of the Maestrello and 
McDaid subsonic jet (ref. 52). The potential core is longer than indicated by the experi- 
mental data but the decay rate of the center-line velocity in regions II and III is well pre- 
dicted. Test case 6 is an example of a jet for which 1/ 6 must have a value near 0.07 in 
the shear-layer region to correctly predict the core length. 

Test case 7 .- For the Eggers Mach 2.22 jet (ref. 49), the potential core length is 
well predicted, but the agreement with the center-line velocity decay curve using (f> = 8° 
is not as good as the corresponding comparison in the previous subsonic jet flow in test 
case 6. The prediction for test case 7 is shown in figure 11(a) with a <£ = 89° calcula- 
tion which corresponds to an immediate jump in length scale at the end of the core to the 
region III value. This represents the usual method of applying the mixing length model to 
jet flows (e.g., Harsha, ref. 19). Obviously the prediction of the center-line velocity 
decay as the far field is approached is significantly improved by using a "length scale 
equation” between the near field and far field. A prediction with <p = 0° is also shown. 
Velocity profiles at three locations are shown in figures 11(b) and 11(c). These results 
indicate that the correct profile is obtained when the center-line value is correctly 
predicted. 

Test case 8.- A value of 1/6 of 0.05 and mixing angle of 8° were used to predict 
the high-temperature subsonic data of Heck (ref. 53). The center-line distributions of 
velocity and temperature for test case 8 are given in figures 12(a) and 12(b), respectively. 
As shown in these figures, the velocity potential core length is again incorrect although 
the velocity decay rate outside the core is well matched. A turbulent Prandtl number of 
0.7 predicts the temperature decay rate correctly. 

Test case 18.- The fully developed axisymmetric jet data of Wygnanski and Fiedler 
(ref. 28) provided a test of the far field length scale model in a subsonic flow (test case 18). 
Starting with the given similarity profile at x/D = 60, a location in the self-preserving 
region, the solution was obtained with 1/6 = 0.05. With an initial velocity ratio u^ju e 
of 70, the computed spreading rate of the jet was found to be linear and in excellent agree- 
ment with the data. Figure 13(a) shows the center-line velocity prediction over the range 
of the data. The similarity of the data was retained in the calculations as shown in fig- 
ure 13(b) where x v0 indicates the streamwise distance from the virtual origin of the jet. 
Equation (14) was used to define 6. 

Test case 19.- Test case 19 is a supersonic high-temperature jet (ref. 53) similar 
to the subsonic jet of test case 8. The turbulent Prandtl number was assumed to be 0.7 
and 1/6 = 0.05 and c p = 8°. Figure 14(a) shows that the predicted potential core is too 
long; however, the correct velocity decay rate is predicted as the far field is approached. 
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Figure 14(b) indicates that similar conclusions apply to the center-line temperature dis- 
tribution. If the core length problem could be remedied for compressible shear-layer 
mixing, the center-line variation would then probably be correctly predicted for the entire 
flow field. 

In all the compressible jet cases it should be noted that the near field (where the 
present calculations give least accurate agreement with data) is the only region where 
large density gradients occur. It is speculated that the mechanism responsible for the 
disagreement in this region is the same mechanism operative in the supersonic free shear 
layers. In the present paper it is suggested that this mechanism is related to turbulence- 
induced static-pressure gradients. Other possible mechanisms operable in this region 
are the "low Reynolds number effect" discussed by Birch and Eggers in paper 2 and the 
1/6 adjustment from two-dimensional to axisymmetric. 

Plane Jet in a Moving Stream 

Test case 13 .- The length scale model shown in figure 2 was used to predict subsonic 
test case 13. No initial profile data were available other than the experimentally deter- 
mined momentum thickness and uj^Ug. A step velocity profile was selected so that 
u j/ u e was matched. Value of 1/6 of 0.07 and a mixing angle of 8 ° gave an excellent 
prediction of the center-line velocity distribution as shown in figure 15. Since the exter- 
nal stream had a constant nonzero velocity, the experimental velocity ratio n^ju e could 
be accurately represented numerically, unlike the previously discussed jets issuing into 
still air. 


Coaxial Jet Flows 

The length scale model shown in figure 4 was used for all flows of this type. A 
value of 1/6 of 0.05 was used in regions I and II. All solutions were started at the noz- 
zle exit with the initial boundary -layer profiles. For the AEDC experimental data (test 
cases 10 , 20 , and 21 ) these profiles were approximated by using the estimates of boundary- 
layer thicknesses given in reference 54. The mixing angles $ 2 , and ^3 were 

taken to be equal in all cases. 

Test case 9 .- Test case 9 represents the data of For stall and Shapiro (ref. 55). The 
suggested initial profile, which was adjusted to account for probable pitot probe positioning 
errors, was used to start the solution. The results of the computation using 1/6 = 0.05 
and all mixing angles equal to 8 ° are shown in figure 16. The experimental center-line 
velocity decay in region II is well predicted by using angles of 8 ° although the potential 
core length is unpredicted. However, it should be noted that the core length is a function 
of <p 2 for a given 1/6. As discussed previously, a larger average 1/ 6 is needed in 
the core. 
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Test case 10.- Figure 17 shows the results of applying the coaxial jet length scale 
model to the subsonic hydrogen-air data of Chriss (ref. 56). Both the center-line velocity 
and center-line hydrogen concentration variations for test case 10 are well predicted with 
1/6 = 0.05 and mixing angles of 8°. 

Test case 11.- Figure 18 shows the computational results for test case 11 which 
represents an inner subsonic air jet mixing with an outer Mach 1.30 air jet studied by 
Eggers and Torrence (ref. 57). Mixing angles of 2° gave better predictions of the center- 
line velocity variation than did angles of 8°. However, neither calculation predicted as 
large a drop in center-line velocity as the data indicate and neither calculation correctly 
predicted the rate of velocity increase as the far field is approached. 

Test case 12 .- Good prediction of the center-line velocity decrease for the hydrogen- 
air data of Eggers (ref. 58) was obtained for test case 12 with 1/6 = 0.05 and mixing 
angles of 8° even though this choice of angle resulted in an underprediction of the potential 
core length. The center -line velocity distribution is shown in figure 19(a). The same 
prediction trend is also observed in the center-line hydrogen concentration variation as 
shown in figure 19(b). The turbulent Schmidt and Prandtl numbers were both 0.8. 

Test case 20.- Test case 20 represents the mixing of subsonic coaxial jets as mea- 
sured byChrissaiid Paulk (ref. 54). The best overall prediction with the coaxial length 
scale model was achieved with mixing angles of 6^ as shown in figure 20. The calcula- 
tions with angles of 6° and 8° are shown for comparison; however, as discussed previously, 
1/6 should actually have had a larger average value through the core region. 

Test case 21.- As shown in figure 21(a) the coaxial model again predicted the cor- 
rect center-line velocity decay rate for the subsonic hydrogen-air data of Chriss (ref. 56). 
However, a value of 1/6 of 0.05 gave an overprediction of the potential core length for a 
mixing angle of 8°. An l/6 of 0.07 with all angles (f> ^i.e., 2 ,< t > 2>) e( l ua ^ ® 

gives a better core length prediction. The concentration profile prediction for test case 21 

is given in figure 21(b). 


Wake Flows 

The length scale model shown in figure 3 was used in computing all wake flows. In 
addition, the effect of including the second-derivative term (eq. (5)) was investigated for 
some of these wakes. 

Test case 14.- The predicted center-line velocity distribution is shown in figure 22 
for the two-dimensional wake data (test case 14) of Chevray and Kovasznay (ref. 59). By 
using a value of 1/6 of 0.07 and a mixing angle 0 of 8°, a reasonably accurate overall 

/ 

prediction is obtained. As expected, addition of the second-derivative term |with — = - 
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results in an increase in the mixing rate, the increase giving a more accurate velocity 
rise in the near wake region while overpredicting the data as the far field is approached. 

A calculation with <p = 89°, which corresponds to using no linear variation in length scale 
downstream of the body, is shown to demonstrate the improved prediction in the near wake 
region with the present linear variation. The effect of the angle disappears as the asymp- 
totic region is approached. A calculation with /3 = 4° and no second-derivative term is 
also shown for comparison. 


The deviation between prediction and experiment in the far field is emphasized when 
the data is plotted in terms of l/w^ 2 where w^ is the center-line velocity defect. 


Test case 15 .- Computations of Chevray's axisymmetric wake (ref. 60) with -=0.05 

6 

and j3 = 8° are shown in figure 23. The initial rapid increase of the center-line velocity 
(due in part to the pressure gradient just behind the body) is not predicted even when the 


second-derivative term 


with 



l 

6y 


is included. 


The second-derivative model, how- 


ever, gives better downstream agreement than the basic model. Test case 15 is an exam- 
ple of a flow in which the value of /3 probably should be high initially and lowered to 8° 
a few diameters downstream (as is done automatically by the method of ref. 20), but such 
optimization was beyond the present intent. 


Test case 16 .- The asymptotic region of the plane supersonic wake of Demetriades 
(ref. 61) is well predicted with | = 0.07 and p = 8°. (See fig. 24.) The second deriva- 
tive increases the mixing rate only slightly for test case 16. 


Test case 17 .- The supersonic (M = 3.0) axisymmetric wake data of Demetriades 

(ref. 62) is underpredicted for test case 17 with - = 0.05 and p = 8° even though the 

6 

second-derivative term is included. (See fig. 25.) The calculation was started at — = 17 

D 

a location actually near the end of the transition region; hence these data may be influ- 


enced by the "low Reynolds number effects." Assuming turbulent flow at — = 17 results 

D 

in an underprediction of the data in the turbulent region f— > 2oV 


Test case 24 .- The two-dimensional supersonic wake data of Demetriades (ref. 63) 
include data in the transition region ^602. 9 = 938. 4j. The wake -generating wedge was 

heated so that transition would occur far downstream of the model. The calculation for 
test case 24 was started in the laminar region at — = 183.7. In reference 64 Demetriades 
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gives the experimental streamwise distribution of the intermittency factor on the wake 
center line for the corresponding unheated flow. This distribution was used to compute 
through the given region of transitional flow. Such an approach had previously been 
employed by Harris (ref. 65) who used the streamwise intermittency function of Dhawan 
and Narasimha (ref. 66) in the computation of compressible boundary layers. The results 
shown in figure 26 indicate that the data were slightly underpredicted in the laminar and 
transitional regions with improved agreement in the fully turbulent region. A wake flow 
mixing angle of 8° was used. 


CONCLUDING REMARKS 

Prandtl's mixing length model has been applied in a consistent manner to compute 
the wide range of free turbulent mixing flows selected as test cases in this compilation. 

On the basis of these computations, the following concluding remarks can be made. 

The mixing length constant if 6 was found to be lower in axisymmetric flows than 
in planar flows. With if 6 = 0.05 in axisymmetric flows and 0.07 in two-dimensional 
flows, the calculations compare very favorably with the data for jets downstream of the 
core, wake flows, and low-speed shear layers. 

For the mixing of flows with differing or developing turbulent length scales (such as 
near field to far field transition region in a jet and wake flows developing from boundary 
layers), some method of correctly determining the local length scale should be used. In 
the present paper, a linear algebraic function based on modeling the physical spreading of 
the turbulence was used. Perhaps a more satisfactory approach would be the use of one 
of the various "two -equation" turbulence models of Spalding’s group (paper 11 of this com- 
pilation) where a differential length scale equation is solved. 

The mixing length model was unable to successfully compute the spreading rates of 
free shear layers with large sustained density differences. This inability is tentatively 
ascribed to the assumption of constant static pressure in the transverse direction (the 
conventional quasi-parallel flow assumption). Therefore, a successful calculation method 
for this class of flows (free turbulent mixing flows with sustained large density changes or 
differences) may necessitate including equations for the mean normal momentum and v' 2 . 

As discussed by Birch and Eggers in paper 2 of this compilation, there seems to be 
an important "low Reynolds number effect" in some of the available data. However, fur- 
ther analysis of the data is needed before this effect can be incorporated into prediction 
methods. 
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(b) Variation of 6 with y. 

Figure 4.- Definition of mixing region width for coaxial jet flows. 
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Test Case 1 



constant-density free shear layers. Test case 1. 
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Figure 6.- Effect of Mach number on spreading parameter for two-dimensional 

free shear layers. Test case 2. 




Figure 7.- Prediction of spreading parameter for two-dimensional variable-density free shear layers. Test case 3. 





















Test Case 6 



Figure 10.- Predicted and experimental center-line velocity distribution for a 
subsonic axisymmetric jet issuing into still air. Test case 6. 



Test Case 7 




Test Case 7 
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Figure 11.- Continued. 
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(a) Center-line distribution. 

Figure 13.- Predicted and experimental velocity for the far field of an 
axisymmetric jet issuing into still air. Test case 18. 
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(b) Similarity profiles. 
Figure 13.- Concluded. 
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Figure 15.- Predicted and experimental center-line velocity for two-dimensional 

jet in a moving stream. Test case 13. 




Test Case 9 




Test Case 10 



(a) Velocity. 

Figure 17.- Predicted and experimental center-line distribution for 
subsonic hydrogen-air coaxial flow. Test case 10. 





Test Case 10 





Test Case 12 
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Figure 21.- Predicted and experimental center-line distribution for 
subsonic hydrogen-air coaxial mixing. Test case 21. 









Figure 22.- Predicted and experimental center-line velocity for a 
two-dimensional subsonic wake. Test case 14. 






Test Case 14 



Figure 22.- Concluded. 
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Figure 23.- Predicted and experimental center-line velocity for a subsonic axisymmetric wake. Test case 15. 
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Figure 24.- Predicted and experimental center-line velocity for a two-dimensional supersonic wake. Test case 16. 





Figure 26.- Predicted and experimental center-line velocity for a two-dimensional supersonic wake. Test case 24. 
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Figure 26.- Concluded. 


DISCUSSION 


H. McDonald : Can we regard this paper as a yardstick by which to judge the other pre- 
dictions; in other words, if you cannot do any better than Prandtl's mixing length, then 
why bother? That is the first question I would like for you to try and answer. The sec- 
ond one is, are you advocating a development of Prandtl's mixing length and should this 
work be pursued, or do you think that this is the point at which to abandon the specific 
concept of mean field-direct relation with the turbulence and go on perhaps to the turbu- 
lence kinetic energy equation and multiequation models of turbulence? 

D. H. Rudy: Since a length scale equation has been added to the basic mixing length 
model, the calculations cannot be regarded as baseline calculations in the classical sense. 
The Imperial College paper has that type of mixing length calculation. Our use of the 
mixing length is merely an attempt to better represent the physics of turbulent mixing in 
nonequilibrium flows. The Conference Evaluation Committee is of course best qualified 
to determine what this approach means in terms of other prediction methods. Personally, 
I think that approaches like the various "two-equation" models which incorporate a length 
scale equation appear promising. 

D. M. Bushnell : We use this thing just as a tool to analyze the flows. We've used other 
tools before and the mixing length is a convenient handle to put on the flow to find out if 
it is doing anything surprising. I think it is obvious that you want to continue development 
and each method in the hierarchy has its own application, and I think I'll leave it at that. 

V. W. Goldschmidt : I guess I have a similar question: What is the basis for determining 
whether the model works or not? Is it strictly a comparison with mean velocity or are 
we going to try to look at turbulence stresses as well? I think if we take that daring step 
your conclusions would have to be rephrased a little bit differently. Am I correct? 

D. H. Rudy: Yes, we have examined only mean flow quantities and not the details of the 
turbulence which are of course more sensitive to modeling. Perhaps the Evaluation 
Committee will comment on this but I do agree with your observation. 

H. McDonald: I'm sorry our time doesn’t permit any more questions. I'd like to proceed 
if I may. 


Written Comments 

S. J. Kline: For many decades most workers, including the discussor, have had the idea 
that the mixing length and eddy viscosity were fundamentally different assumptions. How- 
ever, recent developments suggest this is not so; there are two points. First, both are 
closures which depend only on the local mean strain and scaling on the layer width. Sec- 
ond, it follows from this that whenever one is used, there exists an equivalent formulation 
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which is wholly consistent and equivalent within the accuracy of this level of approxima- 
tion. In the case of Rudy and Bushnell's formulation, for example, one can write 


T 



8u 

9y 


9u 

9y 


(i) 


If we also write 


T 



and equate r in equation (1) and equation (2), we have 


( 2 ) 


e = l 2 AU 

0y 6 

and, since l is scaled on 6 (1/5 = Constant), this is equivalent to e = k6 AU. Thus, 
the differences observed in the two cases should arise from the forms and values taken 
for the constants, and not from a fundamental difference between mixing length and eddy 
viscosity formulations. 

P. J. Ortwerth : The measured static-pressure defects which concern the authors have 
been analyzed by myself prior to the formulation of my model. An order of magnitude 
analysis shows that to first order the proper integration of the radial momentum equation 
becomes 


p + (pv) V = p e 

and the normal velocity component remains a second-order term - that is, 


( 1 ) 


pv 2 « (pv) 'v' (2) 

Further, equation (1) has been subsequently verified by Roger Craig in our laboratory by 
comparing pv'2 measured with a hot wire anemometer, Thermo Systems model 1100, 
and measured static pressure on jets. 
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